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A general approach to consider spatially extended stochastic systems with correlations between 
additive and multiplicative noises subject to nonlinear damping is developed. Within modified cu- 
mulant expansion method, we derive an effective Fokker-Planck equation whose stationary solutions 
describe a character of ordered state. We find that fluctuation cross-correlations lead to a symmetry 
breaking of the distribution function even in the case of the zero-dimensional system. In general 
case, continuous, discontinuous and reentrant noise induced phase transitions take place. It is ap- 
peared the cross-correlations play a role of bias field which can induce a chain of phase transitions 
being different in nature. Within mean field approach, we give an intuitive explanation of the sys- 
tem behavior through an effective potential of thermodynamic type. This potential is written in the 
form of an expansion with coefficients defined by temperature, intensity of spatial coupling, auto- 
and cross-correlation times and intensities of both additive and multiplicative noises. 
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I. INTRODUCTION 

A considerable activity concerning the stochastic 
processes addresses the constructive role of fluctu- 
ations of environment being called as noises. An 
incomplete list of such processes includes noise in- 
duced unimodal-bimodal transitions 0, stochastic 
resonance Q, noise induced spatial patterns and 
phase transitions etc. In general case, consid- 
ering a stochastic dynamics, one should deal with 
a problem to account for correlations between ran- 
dom sources. Along this line several special methods 
were developed. Most popular of them are as fol- 
lows: (i) the cumulant expansion method (ii) 
the spectral width expansion method 0, E| ; (iii) the 
unified colored noise approximation where evolution 
equations for both a stochastic variable and a random 
force are combined within unique equation of motion 

nana 

The wide spectra of works are aimed to explore 
an effect of correlations of fluctuations in extended 
systems (see Ref.jj| and references therein). It is ap- 
peared the interplay among noise correlations, non- 
linearity and spatial coupling leads to a special type 
of noise induced phase transition known as reentrant 
transition . Moreover, such reentrance can be ob- 
served even in the limit of weakly correlated noise: if 
the noise self-correlation time r — > 0, the system un- 
dergoes a single reentrant phase transition [l^l ■ In the 
opposite limit of the strong correlated noise (r — > oo) 
a chain of reentrant phase transitions can take place 
A kind of these transitions is continuous or dis- 
continuous in dependence of the self-action part of a 
bare potential. 

Quite peculiar picture is observed in the case of 
several noises whose cross-correlations can arrive at 
the remarkable and counterintuitive phenomena re- 



*Electronic address: alex@ufn.ru 



lated to the problem of reconstruction of phase transi- 
tions. In such a case the stochastic system undergoes 
a chain of different phase transitions with appearance 
of a metastable phase in spite of the fact that bare 
potential does not assume such a phase [l3j. Un- 
fortunately, nowadays we have a scanty collection of 
works devoted to this phenomenon. Perhaps, it can 
be explained by problems in realization of natural or 
computer experiments in extended systems, on the 
one hand, and by absence of theoretical tools and 
methods to perform the corresponding calculations, 
on the other hand. 

Above pointed phenomena force to reconsider be- 
fore developed ideas concerning the theoretical ap- 
proaches of noise induced phase transitions in ex- 
tended systems. Most articles concern the problem of 
above cross-correlations within framework of special 
models which are difficult to generalize for description 
of reconstruction of phase transitions. Therefore, the 
fluctuation induced rebuilding in the system behavior 
is an open question that should be developed to find 
the unique role of stochastic environment influence. 

In this paper, we consider the general situation to 
present the role of cross-correlation contribution of 
two noises into the picture of phase transitions. We 
explore an extended stochastic system which obeys 
the archetypal model of Brownian particle. The ade- 
quate scheme which allows to specify statistical prop- 
erties of the system with nonlinear kinetic coefficient 
in the overdamped limit is introduced. Within a sim- 
plest model with nonlinear damping, drift caused by 
Landau-like potential and two colored (multiplica- 
tive and additive) noises, we show in what a way 
the system can undergo noise induced phase transi- 
tions. We find out phase transitions of both continu- 
ous and discontinuous character are realized as biased 
phase transitions. We obtain phase and bifurcation 
diagrams to elucidate the cross-correlation is one of 
reasons to appear the metastable states inherent in 
first order phase transitions. 

The paper is organized in the following manner. 
Section II is devoted to development of analytical ap- 
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proach to study noise induced phase transitions on 
the basis of kinetic equation for the probability den- 
sity function and cumulant expansion method. In 
Section III, we apply the derived formalism to con- 
sider noise induced phase transitions in the simplest 
case of the Ginzburg-Landau model with both mul- 
tiplicative and additive noises and kinetic coefficient 
dependent of the stochastic variable. On the basis of 
obtained drift and diffusion coefficients we build out 
an effective one-dimensional stochastic process being 
multiplicative but noncorrelated. Related probability 
distribution function combined with self-consistency 
condition allow us to investigate the corresponding 
phase diagram and stationary behavior of the or- 
der parameter (Section IV). Discussion in Section V 
is based on representation of the effective stochas- 
tic process within framework of the mean field ap- 
proach. Such consideration allows us to study noise 
induced phase transitions by analogy with the stan- 
dard Landau scheme. It is appeared the noise cross- 
correlations and strengthening dispersion of damping 
coefficient arrive at transformation of the continu- 
ous phase transition into discontinuous one. Finally, 
main results and perspectives of the work are col- 
lected in Conclusion (Section VI). 



II. MAIN RELATIONS 

Consider a Brownian particle under influence of an 
effective potential ^[^(r)] and a damping character- 
ized by the viscosity coefficient 7(2:). The generalized 
evolution equations for the scalar variable x and the 
conjugate momentum p read: 



mx = p, 
p + -f(x)p : 



ST 



8x(r, t) 



(1) 

+ g^x)Ur,t) (2) 



where m is the effective particle mass, dot stands for 
the derivative with respect to the time t, r is the 
space coordinate; the effective potential is reduced to 
the Ginzburg-Landau form 



T = 



V (x) + ^\Vx\ 2 



dr 



(3) 



with Vq(x) and D > being a specific thermo- 
dynamic potential and an inhomogeneity constant, 
V = d/dr. The last term in Eq.©, where index 
/1 = 1,2 numerates different noises to be summa- 
rized in accordance with the Einstein rule, represents 
Langevin forces which act with amplitudes g^ (x) and 
stochastically alternating functions Neglecting a 
space correlation, we focus on time correlations be- 
tween forces i.e., 

(C M M)C,(r',t')> =5{v-v')C^{t-t'). (4) 

Inserting Eq.J5J) into the result of differentiating 
Eq.Q over time, it is convenient to represent evo- 



lution equation for the quantity x in the form 
mx + j(x)x = f(x) + DAx + g^x)^ 



(5) 



where f{x) = — AVq/Ax is a deterministic force. 

To study statistical properties of the system one 
needs to find the probability density function P = 
P(p, x, t) of the system states distribution in the 
phase state {x,p}. With this aim, we represent the 
system on the regular d-dimension lattice with mesh 
size A£ = 1 and -f t = 7(2^), fi = f(xi), g Mi = g^Xi)- 
Then, the differential equation with partial deriva- 
tives (0 is reduced to the usual differential equation 



Xi ~\- 'JiXi — fi~\~ ~ , ^ DijXj ~\- g^iCn 



2d 



(G) 



due to representation of the Laplacian operator on a 
grid as follows 



2dSi 



(7) 



where nn(i) notices a set of the nearest neighbors of 
the site i. 

By definition, the probability density function is 
given by averaging over noises of the density function 
p(xi,pi, t) of the microscopic states distribution in the 
phase space: 



P(x l ,p i ,t) = (p(xi,Pi,t)). 



To construct an equation for the macroscopic density 
function P — P(xi,pi,t) we exploit the conventional 
device to proceed from the continuity equation for 
the microscopic one p = p(xi,pi,t): 



dp 
at 



0. 



(9) 



Inserting the time derivative of the momentum p 
mx from Eq.© into Eq.©, we obtain 



dp 
dt 



= (C + Af^Ap 



(10) 



where the operators C and are defined as follows: 

£ = Pi - 

m dxi 



_ 7i 

'i Pi I ? 

m 



Cr - d 
» = ^ l dp-- 



(11) 



(12) 



Within the interaction representation, the microstate 
density function reads as 



P 



e~ ct p 



(13) 
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to reduce Eg . (|10fl to the form 



(14) 
(15) 



A standard and effective device to solve such a type 
of stochastic equation is the well-known cumulant ex- 
pansion method Q. Neglecting terms of the order 
0(1Z^), we arrive at the kinetic equation of the form 



L 

j2 f {%mu{t'))&t' 



(p)(*). 



(16) 



Within the original representation, the equation for the probability density (JSJ) reads 

d 



-P(t) = <£+ C^(r) ( e £r A/>^)] At \ P(t) 



(17) 



r 



If the physical time is much more than a correlation 
scale (t ^> Tfj), we can replace the upper limit of the 
integration by t = oo. Then, expanding exponents, 
we derive to the perturbation expansion 



where collision operator 

oo 

c = c (n) > <? (n) = mw (a^™ 

n=0 

is determined through the commutators 

and moments of correlation function 
1 



r n C^(r)dr. 



(18) 



(19) 



(20) 



(21) 



To perform next calculations we shall restrict our- 
selves considering overdamped systems where the 
variation scales t s , x s , v s , j s , f s , and g s of the time t, 
the quantity x, the velocity v = p/m, the damping co- 
efficient 7(2;), the force f(x) and the noise amplitudes 
<7 M (x), respectively, obey the following conditions: 



v s t s 

fs^s 



e 1 > 1, 7s i s 

g s t s 



^> 1, 



(22) 



e- 1 » 1 



e" 1 » 1. 



These conditions means a hierarchy of the damping 
and the deterministic/stochastic forces to be charac- 
terized by relations 



fs/m 



e < 1, 



9s /m 



e<l. 



(23) 



As a result, the dimensionless system of equations 
Q , (0) takes the form 



at 

-=-e 7 (*)« 



(24) 



e- 1 [f{x)+DAx + g^x)^{v,t)] 



Respectively, the Fokker-Planck equation ijTS|l reads 



(25) 



where the operator 



C = e- 1 C 1 + e- 2 C 2 



(26) 



has the components 



dxi 




(27) 



£2 = 7» "5 — Vi 
dvi 



The collision operator C is defined through the ex- 
pressions l|19[) — (|21|l and the operator (|12|l with the 
momentum pi being replaced by the velocity Vi. With 
accuracy up to the first order in e <C 1 the expansion 
<|19|) takes the explicit form 
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1 / #5* 



9a;,; 



+ 



i:) 2 



dxidvi 



+ 0(e 2 ). (28) 



r 



To obtain the usual probability function V(xi, t) we 
consider velocity moments of the initial distribution 
function P(xi,Vi,t) in the standard form [6j 



V n (xi,t)= v?P(xi,Vi,t)dvi 



(29) 



where integration over all set {vi} is performed. 
Then, if we multiply the Fokker-Planck equation 1(25() 
by the factor v n and integrate over velocities, we de- 
rive the following recurrent relations: 



~ -57- + njiV n + e 
at 



OVn+l 

dx. 



D 



f* + 2dE^ I Vn - x 



= n(n - 1) (A/W - 7i M$) - enM^J 



OVn-l 



dg». 
dx. 



+ 0(e 2 ). 



(30) 



At n = 0, we obtain the equation for the distribu- 
tion function V = Po( x ii t) : 



dV 
dt 



, x dVi 
dx,, ' 



(31) 



The expression for the first moment V\ follows from 
Eq. (|30|l where n = 1 and only terms of the first order 
in e are kept: 



2d 



(32) 



dV 2 

dxi 



M, 



(i) 



(I V 



dV 

dxi 



g»i 



'dg v 



dx,. 



The second moment V 2 can be obtained if one puts 
in Ea. H3U|) n = 2 and takes into account only zeroth 
terms of smallness over the parameter e <C 1 : 



Vn 



M, 



(0) 

//// 



7i 



9*9*P- (33) 



At last, the Fokker-Planck equation takes the 
Kramcrs-Moyal form 



dV 
~dt 



_d_ 

dxi 



(PxV) 



d 2 
dx 



2 (piP) ■ 



(34) 



where effective drift and diffusion coefficients are as 
follows: 



(35) 



5pi ^ dx, 



M (0) 
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2 g^igvi- 



(36) 



To proceed the consideration we need to pass from 
the grid representation to a continuous one. In so 
doing, we use the mean field approximation to re- 
place the second term of effective interaction force in 
Eq.®: 



D 



-^^DijXj -► D(rj-x) 



(37) 



where an order parameter 77 = (x) is defined through 
the self-consistency equation 



xP^^dx, 



(38) 



V n (x) is a solution of the Fokker-Planck equation 
1)34(1. Under stationary condition, the relevant dis- 
tribution function has the form 



V v (x) 



z- 1 



exp 



V 2 {x) r I y £> 2 (x') 

\— 00 

where the partition function 



dx' 



(39) 



2,= 



dx 



V 2 (x) 



exp 



V 2 {x') 



dx' 



(40) 



takes care of normalization condition. The equation 
1)38(1 has solutions within the domain bounded by the 
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Newton-Raphson condition 

oo 

/d 



da; = 1 



(41) 



T)=0 



obtained by differentiating Eq. i|38|) over the order pa- 
rameter rj. 



III. MODEL OF CORRELATION BETWEEN 
ADDITIVE AND MULTIPLICATIVE NOISES 

To apply the general results obtained in Section 
ITT1 we consider in details the simplest model of two 
correlated noises being additive and multiplicative in 
nature. Relevant amplitudes are defined as follows: 



la(x) = 1, g m (x) = sign(x)|a;| a 



(42) 



where the exponent is defined as a S [0, 1] and the 
sign function is introduced to take into account a di- 
rection of the Langevin force. We will focus on the 
prototype system concerning the Ginzburg-Landau 
model with the potential 



V (x) 



1 A 

r 



(43) 



where e is a parameter being dimensionless tempera- 
ture counted off a critical value in negative direction. 



In correspondence with the line of consideration \U 
we take up the viscosity coefficient in the form 



II 



(44) 



where positive index a stands to measure the damp- 
ing weakening near bare state x = 1. 

Next, we suppose the noises are to be Gaussian dis- 
tributed, white in space with zero mean and colored 
in time according to the correlation matrix 



/ Zl e -\i-\/Ta £o£m e -|r|/r c 



C(t) = 



(45) 



(7 m 
\ To 



-MAe < p -M/t 



where a a and a m are amplitudes of additive and mul- 
tiplicative noises, respectively, r a and r m are corre- 
sponding autocorrelation times, r c is a time of the 
cross-correlation between the noises. Moments <|21[) 
of both zero and first orders of correlation matrix l|45|) 
are as follows: 



M(°) = 
M« = 



.2 

a 



CT a <J a Ora 
T 2 



7"c(o- m (T Q ) 



T c {a a a m ) 



(46) 



Then, the expressions J3EJ), l|36(l for effective drift and 
diffusion coefficients take the form 



V 1 = \x 2 ~ l\ a { [D(r) -x) + x(e - x 2 )} + aa m \x\ a 1 [<r a T c + (j m T m sign(a;)|a;| a ]} 
+ 2a x(x 2 - I) 2 "" 1 [a a + ( x m sign( a; )| a ;| a ] 2 , 



(47) 



V 2 = {x 2 - l) 2a [a a + ^sign^)^! ] 2 . (48) 

It is worthwhile to notice at 7(2;) = const (a = 0) the 
additive noise can not give a contribution to the drift 
coefficient (|47ll related to the first of moments (|46ll . 



IV. NOISE CORRELATION INDUCED 
TRANSITIONS 

Addressing to the influence of noise cross- 
correlations on the system behavior we start with 
self-consistency equation lp?%|l where the stationary 
distribution is given through the drift and dif- 
fusion coefficients l|47|l . B48JI . It is well known, at 
phase transitions, the symmetry breaking causes the 
ordered state corresponding to the solution 77 7^ 
within the interval of available values of stochastic 
quantity x, while the disordered phase is related to 
i] = 0. In the absence of the multiplicative noise, the 



only reason to break the symmetry of the stochas- 
tic distribution is the interaction force (|37|l which 
plays the role of a conjugate field related to the or- 
der parameter rj. The principle feature of far-off- 
equilibrium systems with colored noise is that the 
symmetry can be restored due to combined effect of 
both the multiplicative noise and the system nonlin- 
earity 0, 0]. Therefore, the reentrant phase tran- 
sition in such systems is appeared. Our aim is to 
demonstrate the course of the phase transition can 
be crucially changed by means of cross-correlations. 

First, we consider the solution of Ea. (|^8fl at dif- 
ferent values of the noise cross-correlation scale t c . 
As shown in Figure ^ in the absence of both noises 
(a a = a m = 0) and coupling (D = 0) the system 
behaves itself in a usual manner being inherent in 
the square-root law (dashed curve with both verti- 
cal derivative in the point of origin and symmetry 
with respect to the e-axis). Such behavior means the 
maxima appearance of the distribution Ij39(l in points 
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±i/e that can be interpreted as standard noise in- 
duced transition of the second order with mean value 
77 = 15]. With switching on noises and coupling 
the situation is changed principally. First, the above 
symmetry is broken to survive the only negative value 
of the order parameter in the limit of small cross- 
correlations (curve 1). Combined effect of correlated 
noises, system nonlinearity and spatial coupling ar- 
rives at the change of the order parameter sign at 
small values r\. Indeed, as seen from curves 2, 3, an 
increase in the cross-correlation time r c shifts weakly 
negative solution into the positive domain causing 
the reorientation transition at the driving parame- 
ter e r . In addition to this transition, an increase 
in the cross-correlation scale arrives at positive so- 
lutions appeared according to discontinuous phase 
transitions that has reentrant nature within domain 
bounded by both lower e c and upper e c boundaries 
(see elliptic form parts of curves 2, 3 where solid and 
dotted lines relate to (meta)stable 0] and unsta- 
ble solutions). With subsequent growth of the cross- 
correlation time above the value T cr related to thin 
solid curve in Figure ^ back bifurcation happens and 
temperature dependence of the order parameter takes 
one-connected character. This means an increase in 
the control parameter e arrives at the (meta) stable 
branch of positive magnitudes 7/ initially (solid curve), 
then the unstable branch (dotted curve) follows from 
the point e c down to e c and finally the negative 
(meta)stable state is merged. It results in formation 
of a hysteresis loop in 77(e) dependence where both 
(meta)stable and unstable states exist to be solutions 
of Eq.(|38J) (see curves 4, 5). Thus, one can conclude 
both reentrance and reorientation phase transitions 
are inherent in systems with colored multiplicative 
noise. 

In FigureElwe plot a phase diagram in (e, r c ) plane 
to show the influence of the noise cross-correlation 
scale on the bifurcation magnitudes of the control pa- 
rameter. It is seen, at small cross-correlation times 
r c , the negative values of the order parameter r) is in- 
herent in the whole e-domain denoted as N in Figure 
13 With increase in r c the reorientation phase tran- 
sition occurs into the state P related to the positive 
solutions rj > 0. The line of this transition is deter- 
mined by the self-consistency equation (|38|l at con- 
dition i] — 0. At a magnitude To, a doubly bounded 
domain R appears to relate to the reentrant transi- 
tion. At the critical value T cr corresponding to the 
bifurcation curve in Figure^ the domain R passes to 
a region M where all stable, metastable and unsta- 
ble phases take place. Over the critical correlation 
time (r c > r cr ), the dotted curve relates to the lower 
critical value e c of the control parameter in Figure ^ 
(curve 4). 

The influence of the multiplicative noise intensity 
af n on the phase transition picture is demonstrated 
in Figure [3] The upper panel relates to moderate 
cross-correlation times r c where an increase in er^ 
transforms two-connected 77(e) dependence into one- 



connected one varying more fast. With growth of the 
time scale t c , the increase of the multiplicative noise 
intensity arrives at the shrinking the metastability 
domain (see lower panel in Figure |3J). Thus, we get 
the conclusion about dual role of the multiplicative 
noise: at small intensities cr^ <C 1, main influence is 
rendered by the cross-correlations between additive 
and multiplicative components of the noise to sharpen 
the phase transition (see Figure^); on the other hand, 
a raising the intensity of the latter component up to 
of n ~ 1 smears this transition. 

According to Figure 0] only one difference of the 
phase diagram in the axes (e, D) from the situation 
depicted in Figure takes place. Here, at small mag- 
nitudes of the control parameter e, the domain P of 
the positive valued order parameter 77 > relates to 
the whole region of the coupling parameter D — con- 
trary to low bounded domain of the cross-correlation 
times t c . 

To find relations between noise exponent a and the 
control parameter e we consider the phase diagram 
in (e, a) plane. It is appeared for noises with weak 
cross-correlation (r c — > 0) the new phase arises only 
at small enough values of a which define the power 
of the multiplicative noise (see dashed curve in Fig- 
ure EJ. In other words, considering the class of sys- 
tems with both additive and multiplicative noises, 
one should mean that ordering processes are possi- 
ble in the case of weak cross-correlation only if the 
multiplicative noise has a weak power. For the sys- 
tems with a — > 1 weak cross-correlation can not in- 
duce new phase formation. According to solid-dotted 
curves in Figure an increase in r c leads to appear- 
ance of the e small valued domain where the reorien- 
tation transition takes place with a-growth. Besides, 
domains of both positive and negative order param- 
eters, being reoriented, join with metastable phase 
region at small values of index a. However, an in- 
crease in a leads to the reentrant phase transition for 
long range cross-correlations. At small and moderate 
values of the noise exponent a the system behavior is 
inherent in the hysteresis loop formation. 

The influence of the damping exponent a ^ on 
the breaking symmetry picture is shown in Figure El 
It is appeared an increase in a transforms the 77(e) 
dependence in a manner similar to the influence of 
the cross-correlation time r c . Indeed, passage from 
two-connected 77(e) dependence related to curves 1 to 
one-connected curves 2, 3 can be provided with both 
r c increase and a growth — quite similarly to the 
77(e) dependencies variation in Figure ^ The conclu- 
sion about similarity of the influences of the damping 
exponent a ^ and the cross-correlation time r c is 
confirmed with phase diagram in plane (e, a) which 
is topologically identical to the same in Figure at 
small correlation times (see Figure Et) - According 
to the Figure [7jp the cross-correlation time increase 
arrives at the reentrant phase transition (within the 
domain R) due to appearance of additional region N 
related to negative values of the order parameter. 
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Finally, we set up the properties of the Langevin 
sources which allows us to produce the ordering pro- 
cesses in the system. With this aim, we plot the cor- 
responding phase diagram in (r c , a) plane (see Figure 
IHJ). Here, the system undergoes reorientation transi- 
tion related to the transforming the negative valued 
order parameter into the positive one if t c increases. 
On the other hand, at small e, increasing the expo- 
nent a of multiplicative noise, we can make the system 
to undergo a chain of phase transitions at which the 
parameter 77 changes the sign three times as maxi- 
mum (at s = 6.5 and r c = 2.85 for example, see Fig- 
ure[Bp>)- The physical situation becomes more simple 
with an increase in e (Figure |Sp). 



V. DISCUSSION 

To understand main features of the system under 
consideration we proceed from equation of effective 
motion 



x = Vi{x) + y/T>2(x)Z{t) (49) 

related to the Fokker-Planck equation 1(34(1 . In differ- 
ence of the initial noises £p (t) in Eq. JSJ effective noise 
£(t) is of white-type: (£(t)> = 0, (£(t)£(0)) = S(t). 
Within the mean field approach, Ea. 1(49(1 takes the 
form 

dF 

r ] ={V 1 {x)) — , F = AF{ V ) - h V . (50) 

With accounting definition 147(1 . where the sim- 
plest set of indexes a — 0, a — 1 is chosen, a 
thcrmodynamic-type potential AF(r]) and a field h 
are defined by the following expressions: 

AF{r } ) = - e -±^r? + \r 1 \ e m ^r m a 2 m ; (51) 



h = T c <T a a m . (52) 

Comparison of the first of definitions ((51(1 with the 
bare potential 1(43(1 shows the multiplicative noise ar- 
rives at increase of the control parameter e due to 
addition e m whose magnitude is proportional to the 
noise intensity a 2 ^ with the coefficient r m being self- 
correlation time. As a result, a growth of the multi- 
plicative noise intensity in Figure [3] causes increasing 
the order parameter 77 at small magnitudes of the 
control parameter e. A smearing of the related de- 
pendencies 77(e) at moderate values e is caused by 
effective field h. 

This field is inherent in the cross-correlation effect 
fixed by the characteristic time r c and intensities a a , 
a m of both additive and multiplicative noises. Ac- 
cording to Ea. 1(50(1 the field h leads to deepening the 
right minimum of the thermodynamic potential F(rj). 



If cross-correlation effects are so slight that the con- 
dition 

r c <C ( £ + £ ") 3/2 ; c ^ 0.385 (53) 

is applied, the field h is less than a critical value 
h c = C(e + e m ) 3 / 2 and the right minimum of the 
thermodynamic potential F(r)) has a local character. 
It means the positive order parameter appears within 
a two-bounded interval e c < e < e c (see curves 2, 3 
in Figure^. With strengthening cross-correlations, 
when the condition h < h c ceases to be valid, a bar- 
rier between right and left minima disappears and 
a domain of definition of the positive order parame- 
ter becomes bounded by the only upper boundary e c 
(curves 4, 5 in Figure H). 

With passage to the general case o / 1, the ther- 
modynamic potential in Eas. ((50|) . ((51|> takes the form 

F = AF( V ) - h signal , (54) 

AF = - (|t7 2 + ) + It; 4 (55) 

that differs from the initial one by replacement 77 — > 
sign(?7)|?7| a . As a < 1, this replacement derives to 
more strong variations of the thermodynamic poten- 
tial V(ji) within the actual domain 77 < 1 that can 
arrive at the appearance of local minimum at mod- 
erate values 77. As a result, a decrease of the index a 
derives to metastable phase — in perfect accordance 
with Figure 

To ascertain the effect of the index a we consider 
mean field approach in the extreme case a = 1, a = 1. 
Here, the thermodynamic potential in the equation of 
motion ((50(1 takes the form 

F = F^Thv+AFin), F $ = l- e -±±?L±ih (56) 

where reference points F^_ correspond to the domains 
77 < —1 and 77 > 1, respectively, the field h is deter- 
mined by Ea. ((52(l and the addition AF(rj) is defined 
by the expansion 

AF ee ^ + f T7 3 + ^ 4 + ^ + (57) 
2 3 4 5 6 

with the following coefficients: 

A= =F(s + £ m ) + 2ct 2 , B ee (A±T c )a a a m , 

C = ±[l + (e + e m )]+2(a 2 m ~a 2 a ), (58) 

E ee -4a a a m , G = =pl - 2of„ ; 

the upper and lower signs ± in the first Ea. 1(56(1 and 
Eas. (|55|) relate to the domains |?7| < 1 and \r)\ > 1, 
respectively. Comparing these equations with the po- 
tential 1(51(1 addressed to the index q = 0wc convince 
the dispersion of the damping coefficient 1(44(1 arrives 
at transformation of the second order phase transi- 
tion into the first one — as it follows from Figures |H| 

m 
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The form of the thermodynamic potential given by 
Eqs.JJ2"J — (|58|l is shown in Figure El as a function of 
the order parameter. It is seen this dependence has 
three well pronounced minima to be inherent in the 
first order phase transition. Moreover, comparison 
of the curves 1 and 2 shows the switching on the 
interaction field (|37fl arrives at the gradient of the 
dependence F(x) inducing the break symmetry. On 
the other hand, if we would like to pass from the Ito 
calculus above used to the Stratonovich one, we have 
to add the term ^T>' 2 (x) to the drift coefficient T>i{x) 
in Eq. I|f)()[) Then, the thermodynamic potential 
F = — jT>i(?])dr] obtains the addition — i2? 2 (?7) to 
transform the potential given by Eqs.lO — C2§ to 
the form 

F = F § - hr) + AF{r,); 

% = g ^ ± ^r c a a a m , 

h = (l± r c )a a a m , 

AF = - v 2 + -77 3 + - v 4 + -rf + -v 6 , (59) 

A = T£ + [^l - (1 ± r m )a 2 m ] , 
B = (10 ±T c )a a a m , 
C=±(l + e)- [4a 2 , - (6 ± r m )a 2 m ] , 
E = -9cr a a m , G=Tl-5cr^ 

where we take into account Ea. H48|) at a = 1, a = 1. 
Comparison of the corresponding dependencies F(r]) 
shown in insertion of Figure [5] with initial ones F(r]) 
depicted in the main panel shows the Stratonovich 
addition promotes to transformation of the first order 
transition into the second one. 

We proceed with consideration of the form of the 
probability distribution function Ij39|l that is respon- 
sible for the reorientation transition related to curve 
1 in Figure |2t>. As shows comparison of the curves 
a and e in the Figure ITUk . positive magnitudes of 
the order parameter r\ > is related to the dis- 
tribution whose right maximum has a larger height 
and is a wider than the left one (and vice versa 
at -q < 0). Much more complicated picture takes 
place with growth of the correlation time r c when 
strongly pronounced maximum of the distribution 
(|39fl is transformed from the left into the right one 
by means of passage via the bimodal dependence (see 
curves f3, 7 and i5 in Figure ITOb ). 



Above considered situation is picked out to address 
to the constant damping coefficient (|44|l when the 
distribution (|39fl has smooth form due to the index 
a = 0. With passage to general case «/0 the depen- 
dence Vri(x) obtains the pair of the strong maxima in 
symmetrical points x = ±1 (see Figure 1111) . The 
analytical form of these maxima follows from the es- 
timations 



T>i(x) ~ 2a (a a ± a m ) x(x 2 - 1) 
V 2 {x)~ K ± <r m ) 2 |z 2 - 1| 



2a-l 



(60) 



that are given by the dependencies l|47|) . (|48|l near 
the points x = ±1. As a result, we arrive at the 
integrable singularities 



P v (x) 



z- 1 

v 



(a a ± fi m ) \x 2 - 1| 



(61) 



which have the form of the maxima shown in Figure 
1111 In contrary, near the point x — one has the 
estimations 

V 1 {x)~aa a a m T c \x\ a ~ 1 , V 2 {x) ~ o 2 a (62) 
which arrive at the expression 

V [x) ^ ^ exp (^ Tc \xA ~ (63) 

Thus, the singularities of the drift coefficient T)\ (x) 
at the point of origin has integrable character to de- 
rive to the finite value of the probability distribution 
function. 

Traditionally, one is taken to present the distribu- 
tion function l|39(l in the Boltzmann-Gibbs exponen- 
tial form 



(64) 



where an effective potential 



V ef = -a 2 a J |i||dx + a\ ln2Mz) (65) 

is introduced to govern by the probability distribution 
in the usual manner . Usage of the definitions l|47ll 
and in the simplest case a = 1, a = derives to 
explicit form of the potential (|65|) : 
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V ef = -Hx + V(x); H = Dr]+(T c -2)a a a m , V = -x 2 + -x 3 + -x\ 
A = (D - e) - (2 + r m )a 2 m + 2—(Dr 1 + r c a a a m ), 

B = 2^ [(e - D) + (1 + r m )a 2 m ] - 3 (^) (Drj + r c a a a m ), 

a \ ' a / 

C = l-(^\ [3(e - D) + (2 + 3T m )a 2 m ] + 4 fe^ + r c( r a a m ) 



where we kept only the terms up to the fourth order 
of the stochastic variable x. The form of the depen- 
dence (|G5[1 is depicted in Figure IT21 for different sets 
of the indexes a and a. Comparison of the respective 
curves confirms the growth of the index a promotes 
to strengthening pair of the strong minima at points 
T) = ±1 of the bias dependence V e f(x). 



VI. CONCLUSIONS 

In this paper we have considered the effect of the 
ordering of stochastic system with two correlated 
noises. In so doing, we have used a model of the 
system with Landau-like potential Vq(x), subject to 
both additive and multiplicative noises with ampli- 
tude of the last in the form of the power-law function 
\x\ a , a £ [0, 1] and affected by x-dependent damp- 
ing with coefficient j(x) = \x 2 — a £ [0,1]. 
Within the framework of both the cumulant expan- 
sion method and mean field theory, the stationary 
picture of the ordered states is investigated in details. 
We have shown the fluctuation cross-correlations can 
lead to the symmetry breaking of the distribution 
function even in the case of the zero-dimensional sys- 
tem. With introducing the spatial coupling, noise 
cross-correlations can induce phase transitions where 
the order parameter rj — (x) varies discontinuously or 
in reentrant manner. We have studied the specified 
interval of magnitudes of system parameters where 
the ordered phase can be formed. With this aim, 
principle phase diagrams are obtained to illustrate 
the role of the multiplicative noise exponent a, spec- 
tral characteristics of fluctuations (auto-correlation 
time r m and cross-correlation time r c ), amplitudes 
of both additive a a and multiplicative a m noises, ex- 
ponent a of the kinetic coefficient j(x), as well as de- 
terministic parameters being the dimensionless tem- 
perature £ and intensity of the spatial coupling D. 

Above studied picture allows us to generalize the 
theory of the phase transitions to the system with set 
of stochastic forces of different nature. Basing on the 
mean field approach we have shown the system can be 
described through the thermodynamic potential F{rf) 
whose construction differs principally from the bare 
potential Vq(x): so, if the latter has the simplest x 4 - 



I 

form, the former is shown to be of the ry 6 -form. Co- 
efficients of related expansion are obtained to define 
terms of the even powers through the dimensionless 
temperature e, intensity of spatial coupling D, auto- 
correlation time T m and intensities of both additive 
a 2 and multiplicative er^ noises; terms of the odd 
powers are defined through the characteristics r c , a a 
and a m of the noise cross-correlations, respectively. 
Thus, we can conclude the phase transition tends to 
transform its character from continuous to discontin- 
uous due to the noise cross-correlation strengthening. 
This trend displays more strongly with growth of the 
index a whose value determines the dispersion of the 
damping coefficient 7(0;). On the other hand, tran- 
sition from the Ito calculus to the Stratonovich one 
promotes to inverse transformation of the discontin- 
uous transition into the continuous one. 

Obtained results can be applied to a consideration 
of the complex systems which are far-off-equilibrium 
and hold several collective degrees of freedom. As 
shows consideration of three-dimensional Lorenz-like 
system with noises being initially additive in nature, 
usage of the saving principle reduces two of these 
noises to multiplicative ones The physical rea- 

son of such a picture is hierarchical subordination of 
different degrees of freedom. According to our pre- 
vious considerations H3]> [13 a typical example of 
such type takes place in solid state physics where a 
reentrant metastable phase can appear if the matrix 
phase relates to random ensemble of defects of dif- 
ferent dimensions subject to the field of plastic flow 
(driven dislocation-vacancy ensemble). Here, in the 
course of plastic flow different defect structures alter- 
nate one another according to picture of the first or- 
der phase transition. Moreover, structural reorienta- 
tion transitions take place where the sign of the order 
parameter is related to the resulting direction of the 
Burgers vectors of dislocation cluster. One more ex- 
ample of above studied behavior gives reentrant glass 
transition in colloid-polymer mixtures |20l |. 

Note finally all presented results have been derived 
for the system with nonconserved order parameter. 
The perspective of further exploration is to investi- 
gate the system with conserved order parameter. 
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FIG. 1: Dependence of the order parameter 77 on the con- 
trol parameter e at a = 0.8, a = 0, a\ = 4.84, — 0.01, 
r m = 0.01, D — 1.0. Curves 1, 2, 3, 4, 5 correspond 
to t c — > 0, r c = 2.5, 3.0, 5.0, 10.0, respectively. Dashed 
curve relates to bare dependence r\ — iy/e, dotted curves 
correspond to unstable solutions. 
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FIG. 2: Phase diagram in (e, r c ) plane at a = 0.8, a = 0, 
= 4.84, ct^ = 0.01, r m = 0.01, D = 1.0. Domains 
denoted as P, N, R and M correspond to positive and 
negative 77 values, reentrant transition and (meta)stable 
phases, respectively. 
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FIG. 3: The order parameter rj vs. the control parameter 
s at a = 0.8, a = 0.0, D = 1.0, a 2 a = 4.84, r m = 0.01 
and cross-correlation times r c = 2.5 (a) and r c = 5.0 (b). 
Curves 1, 2 relates to the multiplicative noise intensities 
Cm = 0.01 and = 0.25. Points from a to e address to 
corresponding curves in Figure ITHI 
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FIG. 4: Phase diagram in (e, D) plane at a — 0.8, a = 0, 
cr^ = 4.84, = 0.01, r m = 0.01, r c = 2.5. 
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FIG. 5: Phase diagram in (e, a) plane at a = 0, a 2 a = 4.84, 
ct^ = 0.01, r m = 0.01, t c = 2.5, .D = 1.0. Dashed curve 
relates to the limit r c — » 0. 
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FIG. 6: The order parameter r) vs. the control parameter 
£ at r m = 0.01, a 2 a = 4.84, = 0.01, D = 1.0, a = 0.8: 
curve 1 — a — 0.2, r c = 2.5; curve 2 — a = 0.2, r c = 5.0; 
curve 3 — a = 0.7, r c = 5.0. 



FIG. 9: The form of the thermodynamic potential F given 
by Eas. l|56|l — Ij58^ (insertion shows the potential l|59|l . 
curves 1 and 2 correspond to D — and D = 1). 
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FIG. 10: Probability distributions addressed to different 
points of the dependence r/(e) related to curve 1 in Figure 
EJa: curves a, e in panel (a) correspond to the values e = 3, 
r\ = 0.56 and e = 18, rj = —4.18, respectively; in panel (b) 
the control parameter is equal e = 12 and curves (3, 7, 8 
correspond to different magnitudes of the order parameter 
f) = 3.15, r) = 1.14 and r] = -3.32. 
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FIG. 11: Probability distribution related at a — 0.5, a — 
0.5, £ = 2. 
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FIG. 12: The form of the effective potential Ijfi5^ : curves 
1, 2 and 3 correspond to the cases a = 1, a = (e = 10, 
r/ = -2.8); a = 0.5, a = 0.52 (e = 10, r\ = -3.0) and 
a = 1, a = 1 (e = 10, 77 = —4.88), respectively. 



